Femtosecond Coherent Control of Spin with Light in (Ga,Mn)As ferromagnets 
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Using density matrix equations of motion, we predict a femtosecond collective spin tilt triggered by 
nonlinear, near-ultraviolet (~3eV), coherent photoexcitation of (Ga,Mn)As ferromagnetic semicon- 
ductors with linearly polarized light. This dynamics results from carrier coherences and nonthermal 
populations excited in the {111} equivalent directions of the Brillouin zone and triggers a subsequent 
uniform precession. We predict nonthermal magnetization control by tuning the laser frequency and 
polarization direction. Our mechanism explains recent ultrafast pump-probe experiments. 
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Long range magnetic order arises from the interactions 
between itinerant and localized spins in a wide variety of 
systems, such as EuO, EuS, chrome spinels, pyrochlore, 
manganese oxides, or (III,Mn)V ferromagnetic semicon- 
ductors [H, 0| • With ferromagnetic semiconductors one 
can envision multifunctional devices combining informa- 
tion processing and storage on a single chip with low 
power consumption. Fast spin manipulation is of great 
importance for such spin-electronic, spin-photonic, mag- 
netic storage, and quantum computation applications. 

One of the challenges facing magnetic devices concerns 
their speed. The magnetic properties of carrier-induced 
ferromagnets respond strongly to carrier density tuning 
via light, electrical gates, or current While magnetic 
field pulses and spin currents can be used to manipu- 
late spin on the many-picosecond time scale, femtosec- 
ond spin manipulation requires the use of laser pulses 
0, S]- In ultrafast pump-probe magneto-optical spec- 
troscopy, the pump optical pulse excites e-h coherences 
and corresponding carrier populations, whose subsequent 
interactions trigger a magnetization dynamics monitored 
as function of time via the Faraday or Kerr rotation 

The physical processes leading to femtosecond magne- 
tization dynamics (femto-magnetism) are under debate. 
Open questions include the possibility of direct photon- 
spin coupling, the distinction of coherent and incoher- 
ent effects, and the exact role of the spin-orbit inter- 
action. Following the pioneering work of Ref.[3l, many 
ultrafast spectroscopy experiments were interpreted in 
terms of a decrease in the magnetization amplitude due 
to transient thermal effects 0, Si ■ Observations of light- 
induced changes in the magnetization orientation were 
also mostly attributed to the temperature elevation, 
which leads to transient changes in the magnetic easy 



axes [qI-HoI. 11 |. Most desirable is nonthermal magnetiza- 
tion control within the femtosecond coherent I4I tempo- 
ral regime, which promises more flexibility limited only 
by the optical pulse duration. Experiments in ferrimag- 
netic garnets were interpreted in terms of an interplay 
between the inverse Faraday effect [l^l and long-lived 
changes in the rnagneto-crystalline anisotropy [4]. In 
(Ga,Mn)As, Ref.|14| reported magnetization precession 
triggered by changes of magnetic anisotropy on a ~100ps 
time scale due to carrier relaxation, while Ref. [l5| demon- 
strated coherent control of the precession. Recently, 
Wang et.al [5| observed two distinct temporal regimes of 
magnetization re-orientation when (Ga,Mn)As is excited 
by 3.1eV photons. Prior to the ~100ps precession, they 
observed a quasi-instantaneous magnetization tilt, which 
was absent for excitation near the GaAs fundamental gap 
(-1.55eV). 

In this letter we calculate the Mn spin dynamics 
triggered by femtosecond coherent photoexcitation of 
(Ga,Mn)As with linearly polarized light. The joint den- 
sity of states for interband transitions has a strong peak 
in the neighborhood of 3eV, due to A3 Ai excita- 
tions along the eight equivalent directions {111} of the 
Brillouin zone (BZ), the A-edge [l3|- Using a full tight- 
binding calculation of the bands, we show that photoex- 
citation at the A-edge is advantageous for non-thermal 
spin manipulation. This is due to the interplay between 
BZ symmetry, spin-orbit and magnetic exchange interac- 
tions, and coherent nonlinear photoexcitation. We pre- 
dict a femtosecond temporal regime dominated by co- 
herent e-h pairs and nonthermal populations. This is 
followed by a second regime determined by the magnetic 
anisotropy of the F-point hole Fermi sea. We demon- 
strate control of the magnetization reorientation via the 
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photoexcitation frequency and polarization direction and 
by changing the initial easy axis. 

We use the Hamiltonian H{t) = Hh + Hexch{t) + HL{t), 
where Hb — Hq + Hso + iJpd describes the bands in the 
absence of photoexcitation [17]. i/o describes the states 
in the presence of the periodic lattice potential, Hso is 
the spin-orbit interaction, represented by on-site spin- 
dependent terms, while Hpd is the mean field interac- 
tion of the hole spin with the ground state Mn spin Sq 
(parallel to the easy axis) [1, • Here we focus on the 
metallic regime (hole densities ^10^°cm~'^), where the 
virtual crystal approximation applies [2] . To describe the 
high momentum photoexcited states, we diagonalized Hb 
using the Slater-Koster sp^s* tight-binding Hamiltonian 
[19I . We thus obtained conduction electron states created 

and valence hole states created 
k is the crystal momentum 
and n labels the different bands. The holes experience 
an effective magnetic field colinear to So, acting only on 
p-orbitals, that lifts the degeneracy of the heavy (hh) and 
light (Ih) hole GaAs bands by the magnetic exchange en- 
ergy Aprf — f3cS, where c is the Mn density and S the Mn 
spin amplitude The interaction (3 is assumed inde- 
pendent of k. In H-VI semiconductors, a direct theory- 
experiment comparison suggested that (3 decreases 
along the A direction, likely due to the k-dependence of 



eters from the experimental measurements ll|, |2J] : 
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the hybridization with Mn ions [21| . Hexch is the Kondo- 



like interaction between the hole spin and the photoex- 
cited deviation, AS, of the Mn spin S from Sq flS*!. The 
carrier coupling to the optical pulse is characterized by 
the Rabi energies c?„i„k(i) — dmnic exp [— /:^/r^], where 
Tp=100fs is the pulse duration : [12l | 



Ell — Kc{S^Sy 4- S^S^ + SyS^) -l- KuzS^ K^SxSy . (2) 

The total energy Eh depends on the direction of the 
magnetization unit vector S [17]. In contrast to mag- 
netic insulators Q, the local Mn moments in (Ga,Mn)As 
are pure 5=5/2 spins with angular momentum L=0, so 
the localized electrons do not contribute to the magnetic 
anisotropy. Kc is the cubic anisotropy constant, Kuz 
is the uniaxial anisotropy constant, due to strain and 
shape anisotropy, and describes an in-plane uniaxial 
anisotropy observed experimentally [3, [S] . Here we ne- 
glect for simplicity the light-induced temperature eleva- 
tion and assume the equilibrium values of the anisotropy 
parameters jTl, 24|. For sufficiently large Kuz >0, So lies 



within the x-y plane. For 0< Ku < Kc, as is the case at 
low temperatures, there are two in-plane metastable easy 
axes, X~ and 1^+ as in Ref.jB], which point at an angle 
(j), sm2(j) — Ku/Kc, with respect to the x-axis [100] [2J]. 

With '^3eV photons, we excite high momentum states 
along the {111} equivalent directions in the BZ, which are 
well separated in energy from the thermally populated 
states [16|. We can then distinguish between thermal 
[25! and photoexcited carrier contributions to the mean 
field equation of motion of S(t)=So+AS(i), 



dtS = -7S X H 
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SxdtS, (3) 



where 7 is the gyromagnetic ratio, Asj^ is the deviation 
(from its thermal value) of the total hole spin 



^Lit) = -J2 dmnkit)eUhl^^ + h. 
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We consider linearly polarized optical field E(t) propa- 
gating along the [GDI] crystallographic axis (z axis). The 
dipole matrix elements /imnk, d„i„k=/imnkE(t), were ex- 
pressed in terms of the tight-binding parameters [3] by 
considering the matrix elements of Vki?6(k) (2^ . 

In addition to the photoexcited states, we consider the 
effect of the thermal holes. Wang et.al. [§] measured an 
upper bound of ^200fs to the hole spin relaxation time 
in InMnAs. The Fermi sea spin relaxes on a time scale of 
lO's of fs due to the interplay between disorder and spin- 
orbit [13] and other Q interactions. For such fast relax- 
ation, we assume to first approximation that the thermal 
hole spin adjusts adiabatically to the instantaneous S(t) 
[3] and describe the effects of the Fermi sea bath via its 
total energy Eh{S) [13] ■ This Fermi sea populates va- 
lence states close to the F point. In view of uncertainties 
such as the population of impurity bands and the origin 
of strain [3, [ij] , we adopt the general form of Eh dictated 
by the symmetry 17[ and extract the anisotropy param- 
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a is the Gilbert damping coefficient [2y]jand H*'* = 
thermal hole anisotropy field }25|. 
We describe the itinerant spin and charge dynamics 
within the mean field approximation by deriving equa- 
tions of motion for the carrier populations and coherences 
[1^ . The nonlinear e-h coherences are given by 

i(9t(/l_knekm) = (Ckm + ^kn " */^2) (^-knGkm) 



^mnk {t) 



1 - (^-kn^-kn) - (eL^km) 



+ E rfmn'k(i)(^Lkn''*-kn) + ^ dm'nk (t) (fik^/ 6km / 

+f3cY,^m ■ st„- {h-^weurn). (5) 



The nonlinear contributions to the above equation in- 
clude Phase Space Filling [13] (second Hne), coupling to 
carrier coherences (third line), and transient changes in 
the hole states due to interactions with the light-induced 
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FIG. 1: (Color online) Femtosecond Mn spin dynamics for 
hole relaxation faster than the optical pulse, (a) and (b): 
ASz/S of the two in-plane magnetic memory states X~ (solid 
line) and K"*" (dotted line) for two frequencies ujp that excite 
different valence bands. Dashed-dotted curve in Fig[TJb): 
Wp=1.6eV. Symbols: non-thermal dynamics (H*''=0). In- 
set: Experimental ABk/^k 0] of X~ and for u}p—S.leV. 
(c): Zero-k magnon oscillations following (a), (d): Popula- 
tions and inter-valence band coherences (inset) corresponding 
to (b). Apd=130meV, Tz^SSfs, i«'c=0.014meV, Ku=Kc/3, 
Ku,=5Kc, a=0.03, E^2 x 10^ V/cm (d ~ 0.6-2 meV). 



AS(t) (fourth line). The hole populations and inter- 
valence band coherences are determined by the equation 

+ ^ '^*„„k(^)(^-kn'ekm) " ^ dm„'k (t) (/l-knekm) * 



I 

-/3c^AS(i)-sr„, (^^-k;^-kn'), 



(6) 



while similar equations are obtained for the electrons. We 
note that the hole populations are coupled to the inter- 
valence band coherences due to interactions with AS(t). 
rjj„, , n^n\ are the dephasing rates and rj5„=l/Ti the 
population relaxation rate, which describe the photoex- 
cited hole scattering and disorder/defect trapping effects. 

We consider Sg parallel to an in-plane easy axis, X~ or 
Y+ [i] {S^=0 initially [13). Figsl^a) and (b) show the 
femtosecond stage of the Mn spin dynamics triggered by 
a linearly polarized pulse tuned close to the A-edge with 
a fluence ~7/iJ/cm^ as in Ref.0]. An out-of-plane mag- 
netization tilt develops on this time scale, consistent with 
the experimental results for A9k/0k ~ ASz/S (inset), 
where 9k is the Kerr rotation angle. This tilt practically 
disappears for ^1.6eV photoexcitation of correspond- 
ing states close to the Fermi surface (dashed curve in 
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FIG. 2: (Color online) (a): Femtosecond Mn spin dynamics 
for different Ti . (b) Hole population relaxation corresponding 
to (a), (c) and (d): Control, via the optical field polarization, 
of femtosecond dynamics (c) and Zero-momentum magnon 
oscillations (d). d ~ O.lmeV. Other parameters as in Fig[lja). 



Fi g . 1( b)). For the anisotropy parameters of (Ga,Mn)As 
ll|,|2J|, the thermal contribution to Eq.(I21), H*'', plays a 
minor role on the sub-picosecond time scale. When the 
photoexcited hole dephasing and relaxation occurs faster 
than the pulse duration (Figdfd)), AS develops during 
the optical excitation, due to interactions with the co- 
herent e-h pair spin. This initial dynamics is followed 
by a second temporal regime, Fig[TJc), dominated by 
H*'' and characterized by zero-momentum magnon os- 
cillations with frequency 



47^ + Ku 



Kr 
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(7) 



If the hole relaxation is slower than the pulse duration, 
AS(t) develops on a time scale ~ Ti, the population 
relaxation time (Fig[2I^a) and (b)). We conclude that the 
femtosecond magnetization reorientation is governed by 
the dynamics of both coherent and nonthermal holes. 

Due to its nonthermal origin, the femtosecond magne- 
tization tilt can be controlled via the coherent photoexci- 
tation. First, it increases with intensity (compare Figs[T] 
and [2]). The sign of /S.Sz can be controlled by tuning the 
photoexcitation frequency, which controls which bands 
are excited, and by changing the direction of So (com- 
pare FigsUJa) and (b)). Fig [T] demonstrates femtosecond 
resolution of the easy axis direction. Another means of 
control is demonstrated by Fig[21 which shows the depen- 
dence of the tilt (Fig[2fc)) and oscillations (FiglSJ^d)) on 
the direction of the optical field polarization as the latter 
is rotated within the x-y plane. An analogous effect was 
observed in ferriniagnetic garnets [3]. 

Similar to Ref. 18|, our results can be interpreted in 
terms of photoexcitation of a pulsed hole spin compo- 
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FIG. 3: (Color online) Photoexcited total hole spin compo- 
nents parallel and perpendicular to Sq: (a) Full calculation 
for the parameters of Figl2](b) Hso=0. 

nent As'*-'- perpendicular to Sq. Our mechanism should 
be contrasted to the inverse Faraday effect, which pre- 
dicts that non-resonant excitation with circularly polar- 
ized light induces an effective magnetic field parallel to 
the direction of light propagation in a non-absorbing ma- 
terial [l^. The interpetation of Ref. [l3| relies on equi- 
librium concepts such as free energy. Here we develop a 
non-equilibrium theory that treats both interactions and 
coherent nonlinear optical effects similar to the Semicon- 
ductor Bloch equations Our theory describes the 
dynamics for both resonant and nonresonant photoexci- 
tation and takes into account the quantum mechanical 
coherence between all bands. We may interpret AS(t) as 
triggered by a femtosecond effective magnetic field pulse 
that does not point in the direction of light propagation. 

We now turn to the role of the interactions Hpd and 
Hso in photoexciting As''-'-. The hole spin, Eq.(j4]), is 
determined by the density matrices (/iLkn^-kn'), which 
can be simplified by expanding Eqs.([5]) and ([6]) up to sec- 
ond order in the optical field, and by the spin matrix ele- 
ments. Hpd and Hso a-i'e characterized by the magnetic 
exchange, Ap,^, and spin-orbit, Aso '^SSOmeV, energies. 
In samples where Ap,^ ^ A50, the hole eigenstates are 
spin-polarized almost parallel to Sq for all k, i.e. sjj^,j«0. 
At the same time, the coherences {h^_-^^^h^\^n') are sup- 
pressed when £j^„, — far exceeds the pulse frequency 
width, as for large A^,;. Figure [3] compares the calcu- 
lated photoexcited spin components parallel (Asj^) and 

perpendicular (Asjj and As J) to So to the calculation 
with Hso=0- As''-'- is suppressed for Hso—Q and there- 
fore the light-induced femtosecond dynamics is heavily 
influenced by the (Ga,Mn)As valence bandstructure. 

In the opposite limit Aso 3> Apd, the spin of the pho- 
toexcited hole states is almost parallel to k, which points 
along the {111} equivalent directions. The total As'' 
vanishes if all symmetric directions are excited equally. 
However, Hpd and the magnetic anisotropy break this 
symmetry and introduce a preferred direction along Sq. 
The band energies now depend on the projection of k 
on Sq. The linearly polarized optical field introduces an- 
other preferred direction that changes the Rabi energies 
dmnk- As a result, different k states are not photoexcited 



equally. Finally, for sufficiently small Ap^, the photoexci- 
tation of inter-valence band coherences becomes signifi- 
cant. Since in (Ga,Mn)As Ap^ ~ 150 meV is comparable 
to A501 the spin dynamics results from a competition 
between Hpd and Hso that must be treated numerically. 

In summary, we developed a non-equilibrium theory 
of ultrafast magnetization re-orientation in (Ga,Mn)As, 
triggered by the interactions of photoexcited coherences 
and non-thermal itinerant carriers with local Mn spins. 
We predict an initial femtosecond regime of spin dynam- 
ics and demonstrate non-thermal magnetization control. 
Our calculations explain recent experiments [s^ . 
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